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ABSTRACT 

We extend the spectral range of our stellar population synthesis models based on 
the MILES and CaT empirical stellar spectral libraries. For this purpose we combine 
these two libraries with the Indo-U.S. to construct composite stellar spectra to feed our 
(— > ' models. The spectral energy distributions (SEDs) computed with these models and 

the originally published models are combined to construct composite SEDs for single- 
age, single-metallicity stellar populations (SSPs) covering the range AA 3465 - 9469 A 
at moderately high, and uniform, resolution (FWHM=2.5lA). The colours derived 
from these SSP SEDs provide good fits to Galactic globular cluster data. We find that 
the colours involving redder filters are very sensitive to the IMF, as well as a number 
of features and molecular bands throughout the spectra. To illustrate the potential 
use of these models we focus on the Nal doublet at 8200 A and with the aid of the 



> 

\ newly synthesized SSP model SEDs we define a new IMF-sensitive index that is based 

on this feature, which overcomes various limitations from previous index definitions 
for low velocity dispersion stellar systems. We propose an index-index diagram based 
on this feature and the neighboring Call triplet at 8600 A, to constrain the IMF if 
the age and [Na/Fe] abundance are known. Finally we also show a survey-oriented 
spectrophotometric application which evidences the accurate flux calibration of these 
models for carrying out reliable spectral fitting techniques. These models are available 
■ through our user-friendly website. 

^ ' Key words: galaxies: abundances - galaxies: elliptical and lenticular, cD - galaxies: 

evolution - galaxies: stellar content - globular clusters: general - stars: fundamental 
parameters 
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1 INTRODUCTION 

A quantitative study of the stellar content of galaxies and 
star clusters requires the use of the so called stellar popu- 
lation synthesis models. The method consists in comparing 
observables like colours, line-strength indices or full spec- 
tral energy d istributions (S EDs) to the predictions of these 
models (e.g.. lTinslevlll980l) . Furthermore these models are 
widely employed in many other type of studies such as, e.g., 
to provide the necessary templates for kinematic or redshift 
measurement or to determine galaxy stellar masses with the 
aid of the mass-to-light ratios (M/L). Most recent models 
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employ at least three main ingredients, which determine the 
quality of the predictions: a prescription for the initial mass 
function (IMF), a set of stellar evolutionary isochrones and 
stellar spectral libraries. Galaxy observables are predicted 
by adding the contributions of all possible stars, in propor- 
tions prescribed by stellar evolution models. The required 
SEDs, colours or absorption line-strengths of all these stars 
are obtained from the stellar libraries. 

Colours, absorption line-strengths and other observ- 
ables from these models have been widely used in the lit- 
erature for interpreting the integrated light of stellar clus- 
ters and galaxies. In the recent years, it has become com- 
mon practice to use model SEDs, mainly for single-burst 
stellar populations (SSPs), at moderately high spectral res- 
olution. These models provide new means for lifting the 



2 Vazdekis et al. 



main degeneracies hampering stellar populations studies, 
such as the on e between the age and the metallicity (e.g., 
IWorthevl[T993 ). For example these SEDs have allowed us 
to define new indices with enhanced abilities to disentan- 
gle the age from the metallicity ( e.g., IVazdekis fc Arimotol 
19991: ICervantes fc Vazdekisl 120091), or eve n the IMF (e.g., 
Schiavon et al.ll200d : I Vazdekis et ai]|2003h . In recent years 



it is becoming increasingly popular to use these model SEDs 
to derive relevant stellar population parameters, includ- 
ing the Star Formation Histor y (SFH), by means of full 
spectrum- fitting approach (e.g., 



iKoleva et al]|2009l ; Ffbieiro et all 



Cid Fernandes et all [20051 ; 
20111 ). which represents an 



alternative to the more standard approach of fitting a se- 
lected number of line-strength indices (see, e.g-. lTrager et all 
1 19981 . and references therein). Furthermore these model 
SEDs have been used in a variety of applications. For exam- 
ple, these SSP SEDs have been shown to improve the anal- 
ysis of gal axy kinematics for both absorption and emis sion 
lines (e.g.. lFalc6n-Barroso et aill2003l ; ISarzi et al.ll2006l ). 

There are models that ar e fed with t h eoreti- 
cal stellar atmospheres ( e.g.. ISchiavon et al.l 2000; 
iGonzalez-Delgado et al.ll2005l ; ICoelho et al.l I2007T ). whereas 
other models employ empirical libraries, mostly in 



the visible (e.g.. Vazdekisl 19991; Bruzual & Chariot 20031; 


Le Borgne et al. 2004; Conrov & Gunnl20ld; Vazdekis et al. 


201C 


;lMaraston & Str6mbackll201ll; IConrov & van Dokkum 


2011 


), but also in other spectral ranges fe.g..|Vazdekis et al. 


2003 


;IConrov & van Dokkumll201ll; iMaraston & Stromback 


2011 


). Models based on empirical libraries are free from 



the uncertainties in the underlying model atmospheric 
calculations and tend to provide good fit s to both, ab- 
sorpt i on line-strengths and spectr a (e.g. IVazdekis et al.l 

(e.g., 



Maraston et al.l 



2009; 



Maraston & Stromback 2011 and photometric data 



Peacock et ail 120111 ). Unlike 



the models based on theoretical libraries, these models are 
generally restricted in their SSP parameter coverage, as 
most of the stars come from the solar neighbourhood, with 
e.g. a limited metallicity and elemental abundance ratio 
coverage. The quality of the resulting model SEDs relies to 
a great extent on the atmospheric parameters (temperature, 
gravity and metallicity) coverage of the library. 

We published mo del SEDs in t h e nea r- IR at 1.5 A 
(FWHM) resol ution ll Vazdekis et al.l 120031 ) and in the 
optical range (IVazdekis et all l201(t) at FWHM=2.5lA 
|Falc6n-Barroso et al.l l201ll ). which employ the CaT and 
MILE S empi rical stellar spectral libr a ries of ICenarro et ail 
l|200ll ) and ISanchez-Blazquez et al l ([2006), respectively. 
These libraries show improved coverage of the stellar param- 
eters and good flux-calibration quality among other char- 
acteristics. Unfortunately there is a spectral gap around 
8000 A, of about 900 A, which is not covered between these 
two model SED predictions. This gap prevents us to derive 
a variety of widely employed colours, from U to I broad 
band filters, and to exploit potentially interesting absorp- 
tion features in that spectral region, such as the Nal at 
8200 A, which has been shown to be usef ul for constrain- 
ing t he dwarf/giant ratios in galaxies (e.g. iFaber fc Frenchl 
Il980h . 

This is the first paper of a series where we present 
the extended spectral coverage of the SEDs predicted 
by our stellar population synthesis models. For this pur- 
pose we make use of the Indo-U.S. stellar spectral library 



l|Valdes et all [2004) to fill-in this gap and to extend blue- 
ward and redward the spectral coverage of our model SEDs. 
These SEDs constitute an extension of our "base models", 
as we combine scaled-solar isochrones with empirical stellar 
spectral libraries, which follow the chemical evolution pat- 
tern of the solar neighbourhood. The models rely as much 
as possible on empirical ingredients, not just on the stellar 
spectra, but also on extensive photometric libraries, which 
are used to determine the transformations from the theoreti- 
cal parameters of the isochrones to observational quantities. 
Our predicted SEDs are therefore self-consistent, and scaled- 
solar for solar metallicity. In the low metallicity regime, how- 
ever, our models combine scaled-solar isochrones with stellar 
spectra that do not show this abundanc e ratio pattern (e.g., 
lEdvardsson et all 1 19931 ; ISchiavonll2007l ). An extensive com- 
parison of the predicted photometric properties with data of 
globular clusters of the MW and M31 and early- type galax- 
ies from the SDSS is shown in a second paper (Ricciardelli 
et al. 2012; hereafter Paper II). 

The paper is organized as follows. In Section 2 we show 
the main ingredients of our models, including the compos- 
ite stellar spectral library that we constructed to feed these 
models. In Section 3 we describe the reliability and quality 
of the resulting models and the behaviour of the MIUSCAT 
SSP SEDs and colours derived from these spectra. Section 4 
illustrates with various examples the potential use of these 
models. In Section 5 we describe the web-based user-friendly 
facilities to download and handle these model SEDs. Finally 
in Section 6 we summarize our results. 



2 MODELS 

The MIUSCAT SSP SEDS presented here represent an ex - 
ten sion of the models pu blished in IVazdekis et al.l (120031 ) 
and IVazdekis et al.l (|2010h to cover the whole spectral range 
AA 3465 - 9469 A. We provide here a brief summary of the 
main ingredients employed by these models and focus on the 
new addition that allows us to achieve such spectral range 
extension, i.e. the MIUSCAT stellar spectral library. 



2.1 Main ingredients 

We employ the solar-scaled theoretical isochrones of 
iGirardi et ail <|2000h . which cover a wide range of ages, well 
sparsed from 0.063 to 17.8 Gyr, and six metallicity bins 
(Z =0.0004, 0.001, 0.004, 0.008, 0.019 and 0.03), where 0.019 
represents the solar value. The isochrones include the latest 
stages of the stellar evolution, including a simple synthetic 
prescription for incorporating the thermally pulsing AGB 
regime to the point of complete envelope ejection. The range 
of initial stellar masses extends from 0.15 to 7Mq. The in- 
put physics of thes e models was updated with respect to 
iBertelli et al.l |l994) with an improved version of the equa - 
tion of state, the opacities of I Alexander fc Ferguson! (|l994l ) 
and a milder convective overshoot scheme. A helium fraction 
was adopted according to the relation: Y ~ 0.23 + 2.25Z. 

We use the theoretical parameters of the isochrones 
(T e ff, log g, [Z/H]) to obtain stellar fluxes on the basis 
of empirical relations between colours and stellar parame- 
ters (temperature, metallicity and gravity), instead of us- 
ing theoretical stellar atmospheres. We mostly use the 
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met al li city-d ep endent empirical relations of lAlonso et al.l 
1 19961 . Il999h : respectively, for dwarfs and giants. Each 
of these libraries are composed of ~ 500 stars and the 
obtained temperature scales are based on the IR-Flux 
method, i.e. only marginally dependent on model atmo- 
sphere s . We use the empirical comp il ation of ILeieune et al.l 
< |l997f ) ; ILeieune, Cuisinier. fe Buserl (|l998T ) (and references 
therein) for the coolest dwarfs ( T e g ^ 4000 K) and gi- 
ants (T e fr ^ 3500 K) for solar metallicity, and also for stars 
with temperatures above ~8000K. For these low tempera- 
tures we use a semi-empirical approach to other metallicities 
o n the basis of these relatio ns and the model atmospheres 
of iBessell et all |l989l . Il99ll ) and the library of iFluks et ail 
j 19941 ). We also e mploy the metal-dependent bolometric cor- 
rections given by lAlonso et all (|l995l . 1 19991 ) for dwarfs and 
giants, respectively. We adopt BCq = —0.12 and a bolomet- 
ric magnitude of 4.70 for the Sun. 

Several I MFs are consid e red: t he two power-law IMFs 
described in IVazdekis et al.l (1 19961 ) (i.e unimodal and bi- 
modal), the two characterized by its slope /i as a free param - 
eter, and the multi-part power-l aw IMFs of iKroupal (|200ll ) 
(i.e. universal and revised). The ISalpeterl l| 19551 ) IMF is ob- 
tained by adopting the unimodal case with slope fi = 1.3. 
An extensive description of these definitions is given in V03 
(Appendix A). We set the lower and upper mass-cutoff of 
the IMF to 0.1 and 100 M , respectively. 



2.2 MIUSCAT stellar spectral library 

To extend the spectral co verage of our MILES-based mod- 
els (|Vazdekis et al.l |2010| ) we employ three empirical stel- 
lar libraries at moderately-hig h spectral resolution: MILES 
( Sanchcz-Blazquez et al. 2006) and the near-IR CaT library 
of ICenarro et al.l (|200lh . both wit h very good relati ve flux- 
calibration, and the Indo-U. S.FI (|Valdes et al.ll2004l ). which 
covers a significantly wider spectral range. This library al- 
lows us to fill-in the gap left between the MILES and CaT li- 
braries and to extend the wavelength coverage slightly blue- 
ward of MILES and redward of the CaT library. 

MILES includes 985 stars covering the range 3525- 
7500 A with a s pectral resolution of 2 . 5A(F WHM), as re- 
cently shown in Falcon-Barroso et alj (|201 lh . The CaT li- 
brary l|Cenarro et alj|200ll ) consists of 706 stars in the range 
8350-9020 A at resolution 1.5 A. The INDO-U.S. library in- 
cludes ~1200 star s with spectra that typical ly cover from 
3465 to 9469 A. In lFalcon-Barroso et all (|201ll ) we also have 
recently shown that the spectral resolution of this library is 
slightly larger (FWHM=1.36 A) than the approximate value 
provided in the original Indo-U.S. paper (FWHM~1.2A). 

The atmospheric parameters of all the Indo-U.S. stars 
in common with MILES and/or CaT libraries have been 
compared in order to calculate the necessary transforma- 
tions for obtaining an homogenized set of stellar parameters 
with MILES. From this comparison we obtained 



Teff M = 1.02Teff /C / - 108.31 
loggM = 1.0467 log g /t / - 0.1077 
[Fe/H] M = [Fe/H] I£7 - 0.02 

1 http://www.noao.edu/cflib/ 



(1) 
(2) 
(3) 



As it can be seen, a linear correction was required to bring 
the published Indo-U.S. temperatures and gravities onto the 
MILES/CaT system, whereas the metallicities were matched 
by applying a small offset. 

We selected the Indo-U.S. stars whose spectra cover the 
desired wavelength range without gaps. We ended up with 
574 stars, from which 142 have been removed due to the pres- 
ence of strong telluric absorption residuals in their spectra. 
The final catalog is composed of 432 stars. Finally we fiagged 
76 of these stars for showing moderate telluric residuals. 

In order to obtain a constant resolution and dis- 
persion value within the desired spectral range (3465 to 
9469 A) for the composite MIUSCAT spectra, we first con- 
volved the CaT and Indo-U.S. spectra with the required 
gaussian kernels to match the lower resolution of MILES 
(FWHM=2.5 A). The next step consisted in resampling the 
smoothed spectra to match the dispersion of MILES (0.9 A). 

Unlike the MILES and CaT libraries the flux-calibration 
of the Indo-U.S. spectra is not very good as a narrow slit was 
employed during the observ ations to achieve a moderately 
high spectral resolution (see lValdes et al' 1 12004 for details). 

We use the Indo-U.S. spectra to fill-in the gap left be- 
tween MILES and CaT spectral ranges, and also to ex- 
tend toward the blue and red the wavelength coverage of 
MILES and CaT libraries, respectively. For each of the se- 
lected Indo-U.S. stars we use the MILES and/or CaT spec- 
tra, when available, and then plug the relevant Indo-U.S. 
spectral ranges to construct the composite MIUSCAT stel- 
lar spectrum. We find 208 stars in common to these three 
libraries. There are 150 Indo-U.S. stars in common with 
MILES and 114 in common with CaT. If a star is not present 
in either MILES or CaT, or even i n these two libraries , 
we use the interpolators d escribed in IVazdekis et~ai1 (|2010T ) 
and IVazdekis et al.l (|2003T l to synthesize the corresponding 
MILES and CaT spectra, respectively. 

To join the stellar spectra we first calibrate the Indo- 
U.S. spectrum in the region of the gap present between 
MILES and CaT. The relative flux scale of the MILES and 
CaT spectra is set by forcing the spectra to have the V — I 
color corresponding to its stellar parameters, following the 
method described in V03 (Section 3.2.1). The V — I color 
is derived from the photometric relations described in Sec- 
tion 12.11 Note that this method ensures a consistent ap- 
proach throughout the various steps of the model computa- 
tions. The Indo-U.S. spectrum in the gap is forced to overlap 
the MILES and CaT spectra at their red and blue edges, 
respectively. This is performed by approximating the con- 
tinuum shape of the original and re-scaled Indo-U.S. spec- 
tra in the gap with a straight line joining the points at AA 
7360-7385 land AA 8390-8415 A(see TableQJ. The Indo-U.S. 
spectrum is then scaled according to the ratio between the 
flux values on these two straight lines on each of the pixels. 
In the Indo-U.S. edges, i.e. bluewards MILES and redwards 
CaT, the INDO-U.S. spectrum is simply re-scaled to match 
the MILES and CaT continua. We scale the relevant Indo- 
U.S. spectral range by matching the flux of the MILES and 
CaT spectra in well selected narrow bands located at nearly 
the edges of spectral ranges of these libraries. These over- 
lapping bands are chosen to be sufficiently wide to achieve 
the statistics. However these bands are sufficiently narrow 
to avoid prominent absorption features and to enable to use 
most of the MILES and CaT spectral ranges. The limiting 
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Table 1. Adopted merging wavelength ranges for matching the 
stellar spectra of the Indo-U.S., MILES and CaT libraries. 



Libraries Overlap wavelength range (A) 



Indo-U.S., MILES 


3540 


- 3564 


MILES,Indo-U.S 


7360 


- 7385 


Indo-U.S, CaT 


8390 


- 8415 


CaT,Indo-U.S 


8926 


- 8950 



IBIUSCAT 

INDO-US 




-o.z 1 1 1 1 — J — * 1 1 1 1 1 1 1 1 1 1 1 

3.4 3.6 3.8 4 4.2 

Figure 1. We plot against logT c ff the residuals obtained by com- 
paring the synthetic (V — I) colour to the expected colour value 
((V - i")phot) for the MIUSCAT stars. The latter values are ob- 
tained by applying the empirical photometric relations described 
in Section |2 . 1 1 according to the stellar parameters. The synthetic 
colours are obtained by convolving the resulting MIUSCAT (black 
open symbols) and the original Indo-U.S. (red solid symbols) stel- 
lar spectra with the filter responses (Johnson-Cousin system). 
Note the large dispersion obtained with the original Indo-U.S. 
spectra due to the limited flux-calibration quality (see the text 
for details). The solid black line represents the median of the 
MIUSCAT residuals for equally spaced logarithmic temperature 
interval. 

wavelengths of the selected overlapping wavelength regions 
are listed in Table [1] 

Finally as the limiting wavelengths of the spectral range 
covered by each stellar spectrum of the Indo-U.S. library 
differs slightly from each other we consider a common range 
for all the stars (AA 3464.9 - 9468.8 A). 

In Fig. [T]we show the residuals for the synthetic V — I 
colours of both, the resulting MIUSCAT and the original 
Indo-U.S. spectra, compared to the expected colour values 
as a function of temperature. The latter are obtained by 
applying the empirical photometric relations described in 
Section 12.11 which are based on extensive photometric li- 
braries, according to the stellar parameters. The synthetic 
colours are obtained by convolving the stellar spectra with 
the filter responses (Johnson-Cousin system). As expected 
the original Indo-U.S. spectra show large residuals due to the 
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Figure 2. The fundamental parameter coverage of the MIUSCAT 
stars is shown for dwarfs (left panel) and giants (right panel). 

limited flux-calibration quality. The composite MIUSCAT 
stellar spectra show good agreement with the expected pho- 
tometric colours. Indeed for the vast majority of the stars 
we obtain residuals smaller than the photometric zero-point 
(i.e. ~0.02mag). 

2.2.1 Stellar atmospheric parameter coverage 

Figure[2] shows the parameter coverage of the MIUSCAT 
stars for dwarfs and giants (separated at \ogg = 3). The 
parameters of the MIUSCAT stars plotted here are on the 
MILES system, which is the result of an extensive compila- 
tion from the literature, homog enized by taking as a refer- 
ence the stars in co mmon with ISoubiran et al.l (119981 ) (see 
ICenarro et al.| [2007). Fig. [2] shows that, at solar metallic- 
ity, all types of stars are well represented. For metallicities 
lower than [Fe/H] ~ —0.8, however, the coverage is rather 
poor and therefore no reliable models can be computed. For 
the metallicity range between solar and [Fe/H] ~ —0.8 the 
coverage of both dwarfs and giants is good, with the notably 
exception of dwarfs with temperatures above ~7000 K. As 
these stars may represent the hotter Main Sequence, includ- 
ing the Turn Off, the resulting model predictions will be of 
lower quality in comparison to those for the older ages. Fi- 
nally the coverage of metal-rich dwarf and giant stars allows 
us to safely compute SSP SEDs for [Z/H] ~ 0.2, includ- 
ing for the young and intermediate age regimes. A quan- 
titative analysis of the quality of the models based on the 
atmospheric parameter coverage of MIUSCAT is provided 
in Section \2. 31 

2.2.2 U magnitude 

Since the bluest wavelength of the MIUSCAT spectra is 
3464.9 A, we cannot accurately measure the flux in the John- 
son U band, as this filter starts at ~3050 A (|Buser fc Kuruczl 
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1 19781 ). A similar situation applies to the SDSS u filter. As 
these filters are commonly used in the literature and the 
MIUSCAT spectral range provides most, but not all, the 
flux in these bands we have estimated the required correc- 
tions to be able to predict these quantities from the resulting 
MIUSCAT SSP model SEDs. To estimate the m issing flux 
bluewards 3464.9 A we employ the IPicklej (|l998l) library of 
stars, which covers the spectral range A A 1150-25000 A at 
low resolution. Note that this library is mainly composed 
by solar metallicity stars. We estimate the missing fluxes by 
convolving the Pickles spectra with the filter responses of 
these bands to obtain both the total flux and the flux in 
the wavelength range covered by MIUSCAT, i.e, redwards 
3464.9 A. 

The correction factors, given by the missing flux divided 
by the total flux, for the Johnson U filter, fu, are shown 
in Fig. [3] for the different spectral types by adopting the 
two, the limiting wavelengths of MILES (upper panel) and 
MIUSCAT (lower panel). Despite the fact that MIUSCAT 
only extends to the blue 75 A the MILES spectral range, the 
missing flux is reduced by about ~10%. This is not surpris- 
ing given the rapid decay of the U filter transmission at these 
wavelengths. As expected hot stars require larger correc- 
tions, particularly those with temperatures above ~10000 K, 
whilst stars with temperatures below ~4000 K show correc- 
tions factors smaller than 10%. The flux corrections are in 
the range 10 - 15% for the stars that mainly contribute to 
the old stellar populations. We obtain separate fits for dwarfs 
and giants in the MIUSCAT spectral range as follows 

/c/=-7.4540 + 4.0006T - 0.5254T 2 , T < 3.955, dwarfs 
/[7 = -6.1348 + 2.8159T - 0.3106T 2 , T > 3.955, dwarfs 

(4) 

/[7=-9.9324 + 5.3528T - 0.7114T 2 , T < 3.895, giants 
/c/=-7.0157 + 3.2215T - 0.3569T 2 , T > 3.895, giants 

where T means logT c ff. The total U magnitude is derived 
by measuring the flux of the MIUSCAT spectrum through 
the U filter redwards 3464.9 Aand correcting it by the miss- 
ing flux. These fractions are taken into account during the 
integration along the isochrone to obtain the correspond- 
ing missing flux corrections for the SSP. It is worth noting 
that these corrections are strictly valid for solar metallicity 
due to the limited coverage of the iPickles! | 19981 ) library in 
this parameter. However the main contribution to this cor- 
rection factor comes from the missing spectral range, which 
represents a small fraction of the total flux in the U band. 

2.3 MIUSCAT SSP SEDs 

We implement the MIUSCAT s tellar library in the models 
as described in IVazdekis et all (|2010T ). The spectra of the 
stars are integrated along the isochrone taking into account 
their number per mass bin according to the adopted IMF. 
For this purpose, each requested stellar spectrum is normal- 
ized to the corresponding flux in the V band, following our 
prescription for converting the theoretical parameters of the 
isochrones to the observational plane, as described in Sec- 
tion 12.11 Note that the latter is performed on the basis of 
empirical photometric stellar libraries, rather than relying 
on theoretical stellar atmospheres, as it is usually done in 
other stellar population synthesis codes. Therefore the SSP 
SED, S\(t, [Z/H]), is computed as follows: 




0.05 U , , , i , , , i , , , i , , , i , , , i , , , i z 

3.4 3.6 3.8 4.0 4.2 4.4 4.6 

log Teff 



Figure 3. Missing flux correcting factor for the Johnson U fil- 
ter versus logT c ff for dwarfs (crosses) and giants (open squares). 
For this purpose we convolved t he U filter res ponse with the low- 
resolution stellar spectra of the Pickles (1998) library. In the up- 
per panel we show the missing flux factor obtained for the MILES 
spectral range, whereas in the lower panel we show the one corre- 
sponding to the MIUSCAT range. Note that the required correc- 
tion for the MIUSCAT range is significantly smaller. The fits for 
dwarfs (solid lines) and giants (dotted lines) are given in Eq. l4l 



prat 

5 A (i, [Z/H]) = / S\(m,t, [Z/B])N(m,t) x 

J mi 

Fv(m,t, [Z/H])dm, (5) 

where S\(m,t, [Z/H]) is the empirical spectrum, normalized 
in the V band, corresponding to a star of mass m and metal- 
licity [Z/H], which is alive at the age assumed for the stel- 
lar population t. N(m, i) is the number of this type of star, 
which depends on the adopted IMF. mi and m t are the stars 
with the smallest and largest stellar masses, respectively, 
which are alive in the SSP. The upper mass limit depends 
on the age of the stellar population. Finally, Fv{m, t, [Z/H]) 
is its flux in t he V band. The absolute flux scaling scheeme 
is described in lFalcon-Barroso et alJ lj201ll ). 

The spectrum of each requested star, with a given set 
of atmospheric parameters, is synthe sized according to the 
inte rpolating algori t hm d escribed in IVazdekis et al.l (2003) 
and IVazdekis et aH (|2010T ) . This code finds the closest stars 
and weight them according to the distance to the requested 
point (6o, log go, [Fe/H]o) in the stellar parameter space and 
the signal-to-noise of their spectra. The method is optimized 
to minimize the errors in case of gaps and asymmetries in 
the distribution of stars around the requested point. For 
a full description of the algorithm we refer the reader to 
IVazdekis etal] l|2003l ) (Appendix B). 

We also take into account during the integration along 
the isochrone the missing fraction of the total flux in the U 
broad band, i.e. bluewards 3464.9 A, for each star, fu i , which 
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is calculated according to Eq. [4] The resulting fraction for 
the SSP, fu S spj i s then used to correct the flux obtained 
by convolving the MIUSCAT SSP SED with the response 
of the U filter. As a test of consistency of these filter mea- 
surements, we co mputed the fugqp fraction s obt ained from 
the SS P SEDs of iBruzual fc Charlotl ([20031 ) and iMarastonl 
(2005) and applied these corrections to the fluxes measured 
in the MIUSCAT SEDs. The magnitudes computed in the 
two ways are fully consistent, with difference of the order 
of 0.01 mag for all the metallicities. This supports our ap- 
proach of applying the corre ction factors d erived from the 
solar metallicity stars of the iPicklesl |l998) library as dis- 
cussed in Section \2. 2. 2 1 

Once an SSP SED based on the MIUSCAT library has 
been computed we combine it with the original MILES and 
CaT SSP SEDs of the same age, metallicity and IMF, to 
construct the final MIUSCAT model SED. For this purpose 
we follow the same approach described for constructing the 
composite MIUSCAT stars (see Section 12.21 to recalibrate 
the models. This approach allows us to keep the MILES 
and CaT SSP SEDs, which are computed on the basis of 
these two libraries, untouched (except the smoothing ap- 
plied to the CaT SED to match the lower resolution of 
MILES, i.e. 2.51 A). Although the SSP SEDs computed on 
the basis of the composite MIUSCAT stars can be safely 
used in the MILES and CaT spectral ranges, our approach 
warranties a significantly higher quality in these two ranges 
as we rely on model computations involving a larger number 
of stars with better atmospheric coverage. Therefore these 
new SEDs, which are computed on the basis of the compos- 
ite MIUSCAT stars, are only used for the wavelength re- 
gions that are not covered by the MILES and CaT spectral 
ranges. In summary, the MILES and CaT spectral ranges 
contained in the final MIU SCAT SSP SEDs are the ones of 
IVazdekis et~ai1 (120031 ) and IVazdekis et~ai] |201Cl ) SEDs, as 
updated in this work. An example of a composite MIUS- 
CAT SSP SED is shown in Fig.0] Table [23] summarizes the 
spectral properties of these models and the SSP parameters 
for which our predictions can be safely used. The models 
assume a total initial mass of 1 Mq. 

For unimodal IMFs with very steep slopes (fi > 1.3) 
the very low MS stars with M < 0.5 Mq might have a no- 
ticeable impact on the resulting mod el SEDs, particularly 
in the redder wavelength ran ges (e.g. I Vazdekis et al.lll996l ; 
IConrov fc van Dokkum|[20Tll) . In Fig. [S] we show the tem- 
peratures and V band magn i tudes for three different sets 
of models: the IGirardi et all (l2000l) that are employed in 
this work, the ones of ICassisi et al.l l|2000l ) that will be im- 
pleme nted in a forthcoming paper, and those of |Pols et al.l 
<| 19951 ) that were emplo yed in IVazdekis etldl (|l996r T models. 
Note that the models o f IGirardi et al l (|2000D are the hottest 
ones. Furthermore the obtained discrepancy in temperatures 
increases with decreasing metallici ty. The difference in tem- 
peratures between the h ottest, i.e. 
the coolest models, i.e. iPols et al. 

lar to those reported in An et all (120091). when c omparing 
the low-mass stellar models of IGirardi et al.l (|2000l ) to those 
implemented in th e Dartmouth stellar evolution database 
jPotter et al.ll2008h . i.e. as large as ^200 K for solar metal- 
licity. Note however that the obtained differences in the mag- 
nitudes of these stars are much less noticeable as is shown 
in the lower panel of Fig. [5] 



Girar di et all (|2000l ) and 




are rather simi- 
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Figure 4. An illustrative example of a MIUSCAT SED corre- 
sponding to an SSP of solar metallicity and 10 Gyr, computed 
with Kroupa Universal IMF. The joined spectra are shown in dif- 
ferent colours: green for MILES, red for CaT and black for the 
Indo-U.S. Note that in the MILES and CaT spectral regions the 
SEDs are identical to the ones published on the basis of these 
two libraries, except the fact that the CaT model SED has been 
smoothed to 2.51 A (FWHM). 



" Girard' 


et al. (2000) 




" _ _ Cossis 


et al. (2000) 




Pols e 


al. (1995) 




















Z=0.004 






Z=0,008 _ 






Z=0.020 




' ' I I 


''''I | i i i 









y/M 



Figure 5. Effective temperatures (upper panel) and V band mag- 
nitudes (lower panel) of very low Main Sequence stars for three 
different set of models and three metallicities, as indicated in the 
insets of the upper panel. 
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Figure 6. Flux ratio between SSP SEDs, at their nominal reso- 
lution, of solar metallicity and 10 Gyr computed with a Unimodal 
IMF with slopes fi = 1.3, 2.3 and 3.3, adopting two dif f erent mod- 
els fo r stars with M < 0.5 M l|Girardi et al.ll2000l: Ippls et al.l 
Il995h . The spectra were normalized to unity at AA 6466-6480 A. 
See the text for details. 



To explore the net effect on the synthesized SSP SEDs 
of adopting different models for these low-mass stars we re - 
placed the stars with M < 0.5M B of iGirardi et alJ (|200d ) 
by those of I Pols et al.l (|l995l ). i.e, the hottest by the coolest 
models, respectively. All the stellar models were converted 
to the observational plane following the same prescriptions 
described in Section T2TI Fig. [5] shows the flux ratio between 
the SSP SEDs of solar metallicity and 10 Gyr with Unimodal 
IMF of slope 1.3, 2.3 and 3.3, computed with these two sets 
of stellar models. For /i = 1.3, i.e. the Salpeter case, we do 
not find any significant effect. For /i = 2.3 the effects are in 
general smaller than 3%, with the largest impact seen for the 
Na feature at ~5800A. Following the NaD Lick index def- 



inition for this feature (IWorthevlll99l V the obtained index 
difference translates to an observing error corresponding to 
a signal-to-noise around 30. Finally for /i = 3.3 the effects in 
some features can be as large as ~10%, and for that reason 
we do not consider safe such model SEDs. 

To estimate the reliability of the newly computed SEDs 
and the SSP parameter regions where these models can be 
considered safe, we use the norma lized quality parameter Q n 
defined in IVazdekis et al.l l|2010h . To compute Q n for each 
SSP we make use of our algorithm for synthesizing a rep- 
resentative stellar spectrum for a given set of atmospheric 
parameters. The larger the total number of stars and the 
shorter the distance of their parameters to the requested 
ones for each star along the isochrone, the larger the quality 
parameter of the synthesized SSP spectrum is. The obtained 
value is normalized with respect to a minimum acceptable 
value, which comes from a poor, but still acceptable, param- 
eter cove rage of the stellar libr ary. We refer the interested 
reader to lVazdekis et all (|2010T ) for a full description of the 
method. Q n allows us to obtain a quality measure of the 
SSP SED, due to the atmospheric parameter coverage of 
the stellar library feeding the models. 
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Figur e 7. The quality parameter, Q n , defined in lVazdekis et aT] 
(2010) as a function of the SSP age (in Gyr) for different mctal- 
licities, indicated by the line types as quoted in the inset. For 
computing the SEDs we use a Kroupa Universal IMF. SSPs with 
Qn values larger than 1 can be safely used (see the text for de- 
tails). 



Figure [7] shows the value of Q n as a function of the SSP 
age (in Gyr) for different metallicities (different line types) 
and adopting a Kroupa Universal IMF. SSP SEDs with Q n 
values above 1 can be considered of sufficient quality, and 
therefore they can be safely used. Note that these values 
only apply to the spectral ranges not covered by MILES 
and CaT SEDs, which in general provide larger values. As 
expected, the higher Q n values are obtained for solar metal- 
licity. Q n reaches a value of ~5 for stellar populations in 
the range 1-10 Gyr. Outside this age range the quality 
decreases as low-mass and hotter MS stars are less numer- 
ous in the MIUSCAT database (see Fig.0). Both, for the 
supersolar and for the low metallicities we obtain lower, 
but acceptable, values. For [Z/H]=-1.31 we obtain Q n val- 
ues that are slightly below 1, but should not be considered 
safe, particularly for detailed spectroscopic studies that em- 
ploy line-strengths or full spectrum-fitting approach. Finally, 
for lower metallicities the MIUSCAT SEDs are not reliable, 
though these models might still be used for some applica- 
tions such as, e.g., broad-band colours (see Section Note 
that in general the Q n values obtained for the MIUSCAT 
SEDs are lo wer than those of MIL ES SEDs (see second panel 
of Fig. 6 in IVazdekis et al1 l2010l). for which we obtain ac- 
ceptable Q„ values for metallicities as low as [Z/H]= -2.3. 
Table[0]summarizes the safe ranges for the MIUSCAT SED 
spectral ranges that are not covered by M ILES and CaT 
SEDs. For the latter we refer the reader to IVazdekis et al.l 
(2010) and the MILES webpage. 
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Table 2. Spectral properties and parameter coverage of the SSP SEDs 



Spectral properties 



Spectral range 
Spectral resolution 
Linear dispersion 
Continuum shape 
Units 



AA 3464.9 - 9468.8 A 
FWHM = 2.51 A, (<j = 64.1km s" 1 at 5000 A, a = 40.1km s" 1 at 
0.9A/pix (54.0km s" 1 at 5000A, 33.7km s _1 at 
Flux-scaled 
F a /L A- 1 Mq 1 , L = 3.826 x lO^crg.s- 1 



SSPs parameter coverage 



IMF type 

IMF slope (for unimodal and bimodal) 

Stellar mass range 

Metallicity 

Age 

SAFE SSP SEDs for non MILES/CaT spectral ranges 



Unimodal, Bimodal, Kroupa universal, Kroupa revised 
0.3, 0.8, 1.0, 1.3, 1.5, 1.8, 2.0, 2.3, 2.8, 3.3 
0.1 - 100 M 
-2.32, -1.71, -1.31, -0.71, -0.41, 0.0, +0.22 
0.06 < t < 18 Gyr 
(-0.71 ^ [Z/H] sC +0.22 & 0.06 < t < 18 Gyr & ^2.3) 



3 BEHAVIOUR OF THE MODELS 

The behaviour of our model SEDs as a function of the SSP 
parameters in t he CaT and MILES wa velen gth ranges have 
been shown in Vazdek is et al. I (|2003T ) and IVazdekis et al.l 
(2010). We refer the reader to these papers for an exten- 
sive analysis. In this section we focus on the new additions 
from expanding the spectral coverage of the model SEDs. 
Among these new features we discuss here our colour mea- 
surements, which are obtained from convolving the MIUS- 
CAT SSP SEDs with the corresponding filter responses. 

Fig. [8] shows t h e broad-band Johnson- Cousins 
l|Buser fc Kuruczl 1 19781 : ICousinsI Il980h colours, U - B, 
B - V, V - R and V - I, derived from the MIUSCAT 
SSP SEDs for different ages and metallicities, adopting a 
Kroupa Universal IMF with th e zero point set with the 
Vega spectrum of iHavesl l| 19851 ). For comparison we also 
show the p hotometric colours of I Vazdekis et al.l (| 19961 ) . as 
updated in Vazde kis et al.l |2010l ). which are calculated on 
the basis of the relations described in Section 12.11 which 
come from extensive empirical photometric libraries. In this 
figure B — V is the only colour that can be fully measured 
within the MILES spectral range, with no need to use the 
MIUSCAT wavelength range extension. 

In general we see that the obtained residuals are as small 
as typical zero-point uncertainties (~0.02 mag) for nearly all 
the metallicities for the B — V , V — R and V — I colours. 
For the B — V colour this can be attributed to the very high 
flux-cali bration quality of the MILES s tellar spectra as was 
shown in lSanchez-Blazquez et all ((2006). On the other hand 
the good agreement obtained for the V — R and V — I colours 
should be attributed to our method for scaling and joining 
the different spectral ranges. Finally the largest residuals are 
obtained for the U — B colour. Note however that the mea- 
surement of the U magnitude is very sensitive to the various 
filter definitions for this band as is shown in Paper II. Our 
photometric U — B pr edictions are mainly obtain ed with the 
aid of the relations of lAlonso et al. I (|l996l . ll999l ). as for the 
other filters. However these U — B relations are not com- 
pletely homogeneous with the relations obtained by these 
authors for the other filters. We refer the interested reader 
to these authors for a detailed description of the methods to 
derive these relations. 

Another interesting aspect that can be studied with the 



newly synthesized MIUSCAT SSP SEDs is the variation of 
the colours and spectral features with the IMF, due to the 
expanded spectral range. The first row of panels of Fig. [9] 
shows the ratio obtained by dividing the SSP SEDs of solar 
metallicity and various IMFs with respect to the SED com- 
puted with a Unimodal IMF of slope 1.3 (Salpeter case), 
for two different ages (2 and 10 Gyr). All the spectra are 
smoothed to 5 A (FWHM), i.e. the constant resolution that 
characterizes the LIS-5.0A s y stem of index measurements 
introduced in lVazdekis et ail (|2010T ) . The first panel shows 
that the models constructed with a Kroupa Universal IMF 
is rather similar to the Salpeter case for the two ages. The 
larger, though small, difference is only seen for redder wave- 
lengths for the 10 Gyr old model. This colour trend is further 
accentuated for the IMF with the largest slope (fi = 2.3), 
i.e., the larger the slope the larger the effect is. Also, the 
effect is more pronounced as the age increases. 

This effect is originated by the large contribution of 
low-mass MS stars in the steep IMF models. Being their 
SEDs red, the resulting total flux of the model is redder than 
models with flatter IMF slopes. Unlike the blue colours, the 
redder ones a re clearly sensitive to t he IMF. A similar result 
was shown in I Vazdekis et al l jl996l ). 

Apart of the colours we also see in the bottom 
panels a significant IMF-sensitivity for various absorp- 
tion spectral features and molecular bands. Among the 
most prominent features we identify the CaH&K around 
~3950A, Cal at ~4227A, the Nal doublets at ~5900A 
and ~8200A, R a and the Call triplet around ~8600A. 
The IMF-sensitivity of all these features, and other fea- 
tures in redder wavelengths, h ave been recently discussed in 
IConrov fc van Dokkuml (|201ll ) (see references therein). We 
also find emphasized IMF-sensitivities in the TiO molec- 
ular bands redwards 5900 A. These include the prominent 
features seen at ~5950A, ~6200A, ~6600A, ~7200A and 
~8900A. The first two TiO bands are included in the Lick 
system of indices (i.e., TiOi and T1O2) (Worthey94). In- 
dices for the other TiO bands have b een discussed and mod - 
elled for stellar popula tions in e.g., ISchiavon et al.l (|l997t ). 
iFaber fc French! (| 19801 s ) and IVazdekis et~ai1 i|2003l ). More- 
over the slope of the SSP SED around the Call triplet at 
^8600 A has also been shown to be sensitive to the IMF 
slope (| Vazdekis et alj|2003l : ICenarro et al.ll2009l ). In the lat- 
ter paper this slope has been modelled by defining a specific 
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Figure 8. We plot the synthetic broad band U-B B-V V - R and V - 1 colours (black lines) derived from the MIUSCAT SSP SEDs fo r 
different ages and metallicities (as indicated within the panels) and the Kroupa U niversal IMF with the z ero point set wit h Ea ves (1985) 
Vega spectrum. The green lines represent the photometric colours computed by IVazdekis et al, 

I dl996h . as updated in IVazdekis et al,l 

(2010), which are calculated on the basis of the empirical photometric relations described in Section |2, II The obtained residuals are also 
shown for each colour. 



index, named sTiO. An index-index diagram, which includes 
this inde x and the Cal l tripl et at ~8600 A, has been pro- 
posed in ICenarro et ail <|2003h to diagnose the IMF slope. 

It is worth noting that whereas the Call triplet de- 
creases with increasing IMF slope, the nearby Nal doublet 
at ~8200 A increases. Therefore the study of these two fea- 
tures constitutes a powerful diagnostic for constraining the 
IMF. The latter feature, which falls in the spectral region not 
covered by the MILES and CaT models, represents a very 
important improvement of the newly synthesized MIUSCAT 
SSP SEDs. We focus on this aspect in Section 14.21 The sen- 



sitivity of the Call triplet to the I MF in SSP models wa s 
extensiv ely shown and discuss ed i n IVazdekis et al.l ll200&t). 
See al so ISchiavon et al.l (|200dh and IConrov fc van Dokkuml 
i|201ll) . The IMF-se nsitivity of the Nal do u blet at ~ 5900 A 
was a lso s hown in Vazdekis et all lll996t) . ISchiavon et al.l 
|200fj|) and lConrov fc van Dokkuml |201ll) . 

As expected, the flux ratio features that are seen at 5 A 
resolution are generally less pronounced when the spectra 
are smoothed to match the lower resolution of the LIS-14.0A 
system, i.e. FWHM=14A, as shown in the second row of 
panels of Fig. [9] This extent shows that the resolution of 
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Figure 9. Flux ratio between SSP models with various IMFs 
(Kroupa Universal in the left panels, Unimodal with slope fi = 0.3 
in the center and Unimodal with slope fi = 2.3 in the right) with 
the corresponding model with Unimodal IMF and fj, = 1.3 (i.e. 
Salpeter). The first row of panels shows models with solar metal- 
licity and two ages (2 and lOGyr), as indicated in the legend of 
the left panel. All the model s were smoothed to match the res- 
olution of the IVazdekis et al. LIS-5.0A line-index system, 
i.e. FWHM=5.0A. The spectra were normalized to unity at AA 
~6466-6480 A. We keep the same scale for all the flux ratios. The 
second row of panels shows the same flux ratios but with the SSP 
SEDs smo othed to match t he sp ectral resolution of the LIS-14.0A 
system of IVazdekis et al.l <201Cft . i.e. FWHM=14A. Finally, the 
third and fourth row of panels show the flux ratios of SSP SEDs 
with lOGyr and two metallicities ([Z/H]=0.0 and [Z/H]=-0.7), 
as indicated in the legend of the first panel, with all the models 
smoothed to match the LIS-5.0A and LIS-14.uA systems, respec- 
tively. 



the data or the velocity dispersion of the galaxies should be 
taken into account when studying the IMF effects. 

In the third and fourth row of panels of Fig. |§]we fix the 
age of the SSP SEDs to 10 Gyr in order to study the IMF 
sensitivity as a function of the metallicity at the resolutions 
of the LIS-5.0A and LIS-14.0A systems, respectively. Overall 
we see that most of the flux ratio features tend to be less 
pronounced in the metal-poor model. Interestingly, this is 
not the case for the Nal doublet at ^8200 A. 

In Fig. [10] we show the difference in B — V and V — 
I colours with respect to the corresponding values for the 
models with Unimodal IMF with fi = 1.3 as a function of 
the IMF slope, for the two ages and two metallicities. As it 
can be inferred from Fig. [9] which shows the larger IMF- 
sensitivity of the redder spectral regions, the V — I colour 
differences are significantly larger than those obtained for 
B — V . We also see that the main increase of the colour 
difference takes place for IMFs with /i > 2. Finally, a similar 
behaviour is found for the colours of the SDSS filters as 
shown in Fig. 1111 



Figure 10. B — V (left panels) and V — I (right panels) colour 
difference with respect to the SSP models with Unimodal IMF 
with /i = 1.3 (i.e. Salpeter), as a function of the IMF slope. In the 
upper panels two ages (2 and 10 Gyr) are considered, whereas in 
the lower panels we show the results for two metallicities ([Z/H]=- 
0.7 and 0). 




Figure 11. Same as in Fig. UOI but for the g — r (left panels) and 
r — i (right panels) colours of the SDSS system. 

4 APPLICATIONS 

In this section we show several applications to illustrate 
the potential use of MIUSCAT SSP SEDs. We only focus 
on those colours and features that we were not able to 
measure with the MILES/CaT SEDs but are now possi- 
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ble with the present extension of the spectral coverage of 
the models. These examples include broad-band colours, a 
newly-defined line-strength spectroscopic absorption index 
and survey-oriented spectrophotometric applications. 



4.1 Galactic globular cluster colours 

Fig- El shows a comparison of the V — R and V — I colours 
with the Milky Way glo bular c l uster data with E(B — V) < 
0.1, from the catalog of lHarrisI (| 19961 ) . The colours are plot- 
ted against metallicities from the same catalog. All the 
models have ~12Gyr and adopt a Salpeter IMF. We pro- 
vide two flavours of m odels: the pho tometric predic t ions o f 
Vazdekis et al. | (|l996h . as updated in I Vazdekis et al.l £2010), 
and the synthetic colours derived from the MIUSCAT SEDs. 
We see that these two sets of predictions agree very well 
each other, matching the observations for the expected ages 
and metallicities according to independent measurements 
from CMD studies. We also show for compari son t he colours 
from the models o f lBruzual fc Chariot! (|2003ft and lMarastonl 
(2005). The latter provide t oo red values in comparison to 
our predictions and the data. lMaraston fc Strombackl (|201ll ) 
have shown that models employing theoretical stellar li- 
braries tend to provide such redder colours in comparison 
to those based on empirical libraries, as in the present work. 

Further colour comparisons for globular clusters in the 
Galaxy and in M 31 are shown in Paper II. Both colours in 
the Johnson-Cousins and in the SDSS systems are discussed. 
We also show in that paper an extensive and comprehensive 
study for the early-type galaxies of the SDSS out to 2 ~ 0.5. 



4.2 Line-strength index definitions 

We illustrate in this section the potential use of the MIUS- 
CAT SSP SEDs for defining new line-strength indices. For 
this purpose we focus on the most interesting feature that 
can be found in the spectral region that was not covered 
by the MILES and CaT SEDs, i.e., the Nal doublet at 
~8200 A (see Section^. This feature, together with the FeH 
Wing- Ford band absorption at 9916 A, has been recently 
used to suggest a dw arf-enriched scenario for massi ve clus- 
ter elliptical galaxies (jvan Dokkum fc Conrovl[2010l ). There 
is also evidence of Nal gradients suggesting concentrations 
of metal-rich dwarfs toward the cen ters of giant galaxies 
|Carter. Visvanathan fc Pickles! fl98rj) . 

We follow here the same scheeme used to define the 
Lick- style indices of the optical spectral range l|Worthevl 
11994 ). We use a feature bandpass and two pseudocontinua 
bandpasses, which are located at either side of the feature, in 
order to establish the flux level of the continuum at the fea- 
ture. These bandpasses are illustrated in Fig. 1131 whereas 
the limiting wavelengths are listed in Table [3] This is a 
modified version of t he index proposed for this feature by 
ISchiavon et al.l (1 19971 ) that avoids the TiO bandhead absorp- 
tion at 8205 A, which becomes increasingly stronger than 
the Nal lines for decreasing tempera tures in giants (see also 
ICarter. Visvanathan fc Pickles! fl986l ). Our index represents 
an i mprovement over previous definitions in this respect 
(e.g.jFaber fc Frenchlll980l.lServen. Worthev fc Briley|[2005l ; 
IConrov fc van Dokkumll201ll ). In fact lAlloin fc Bical (|l989T l 
propose an alternative discussion about dwarf/giants star 




Figure 13. New index definition for the Nal doublet around 
8200 A (NaI8200A). The feature bandpass is indicated in green. 
The continuum flux level (dashed magenta line) is defined by the 
pseudocontinua at either side of the feature (indicated in blue). 
The plotted SED corresponds to a MIUSCAT SSP model of so- 
lar metallicity, lOGyr and Unimodal IMF with slope 1.3 (i.e. 
Salpeter). The spectrum has been smoothed to match the spectral 
resolution of the LIS-5.0A line-index system, i.e., 5 A (FWHM). 



content in the semi-stellar nucleus of M 31, which is based on 
the metallicity, precisely due to the TiO absorption contam- 
ination. Note however that our index definition is certainly 
an advantage when fitting low velocity dispersion galaxies 
as this contamination cannot be avoided in more massive 
systems. Another advantage of the NaI8200A index is that 
we ad opt a significant l y wide r pseudocontinua (~10 A) than 
in the ISchiavon et all l)l997h index (<1 A), which makes it 
more robust against the signal-to-noise of the observational 
data. 

It is worth noting that we do not intend here to pro- 
vide an optimal index definition that requires a specifically 
dedicated paper, which is ou t of the scope of this work. A 
good example can be found in lCervantes fc Vazdekis! |2009), 
where we proposed a new index for the H/3 feature, with an 
unprecedented ability to disentangle the age of the stellar 
populations. Nonetheless, we have tried several index defini- 
tions for this feature, which turned out to be extremely sen- 
sitive to the velocity dispersion smearing of the two absorp- 
tion lines of this doublet. In fact these lines blend in a single 
feature for velocity dispersions around 200 km s _1 . Therefore 
it is important to assess its dependence on this effect if this 
feature were to be used for constraining the IMF. Our exper- 
iments show that the index definition is more robust against 
such smearing if these two lines are included within the fea- 
ture bandpass. Note also that our feature bandpass is well 
isolated from the two selected pseudocontinua bandpasses. 
Fig. 1 141 shows the behaviour of the NaI8200A index against 
the velocity dispersion for models with Unimodal IMF and 
varying slope. The plots show the results for two different 
ages and solar metallicity, but similar patterns are seen for 
other values of these two SSP parameters. As expected, we 
see that the index decreases its strength with decreasing res- 
olution, but this dependence does not vary significantly with 
the IMF slope. 

Although not shown here, we also find that the 
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[Fe/H] [Fe/H] 

Figure 12. Comparis on of V — R (left panel) and V — I (right panel) colours from models of ~12 Gyr and MW globular cluster data 
from the lHarris] d 19961) catalog. We only plot clusters with small extinction values (E(B — V) < 0.1). The colours are plotted against the 
metallicity from the same catalog. The thick magenta line represents an order 2 polynomial fit to the data. All the models are compu ted 
adopting a Salpeter IMF. I n each panel the dotted black line represents the photometric colour prediction of Vazdekis et all lll996h . as 
updated in Vazde kis et al whereas the thick solid g reen line represents t he synthetic colour derived from the MIUSCAT SEDs. 

We also show for comparison the Bru zual fc Charlo ij <200a) and iMarastorj ||2005|) models in long-dashed blue and dot-dashed red lines, 
respectively. 



Table 3. Bandpass limiting wavelengths for the NaI8200A index. 



Feature 


Pseudo-continuum 


bandpasses (A) 


bandpasses (A) 


8180-8200 


8164-8173 




8233-8244 



NaI8200A index is extremely robust against uncertainties 
in the spectrum shape, which might result from, e.g., flux- 
calibration issues. The differences in strengths obtained from 
measuring the index in both an SSP SED with solar metallic- 
ity and 13 Gyr with flux-calibrated response and the same 
model SED after removing its continuum is virtually nil. 
This is expected as the index definition encompasses only 
80 A, i.e. from the bluest wavelength of the blue pseudocon- 
tinuum to the reddest wavelength of the red pseudocontin- 
uum. 

Fig. [15] shows the behaviour of the NaI8200A index as 
a function of age and metallicity for SSP SEDs with Uni- 
modal IMF with two different slopes, fi = 1.3, i.e. Salpeter, 
and steeper, n = 2.3. We use for these measureme nts the 
LIS-14.0A system defined in IVazdekis et"aH (|20ld ). This 
system has a flux-calibrated response curve and a constant 
resolution, FWHM=14A, along the whole spectral range. 
At ~8200 A, where this feature is located, this resolution 
translates to <j~220kms -1 , which is appropriate for mas- 
sive galaxies. Although the NaI8200A index depends on the 
metallicity and age, it is particularly sensitive to the IMF 
slope, as expect ed from Fig. [9] These results are in good 
agreement with ISchiavon et all |200d ). Fig. [15] also shows 



that the sensitivity to the IMF increases with increasing 
ag e and decreas i ng me tallicity. According to the formalism 
of lCardiel et"afl |l998t ) we obtain for this index an error of 
~0.2 with a S/N(A _1 )~30, for an old and solar metallicity 
model. Such error is similar to the index difference obtained 
when varying the IMF slope from 1.3 to 2.3 (see Fig. I15[) . 
This shows that to apply this index for constraining the 
IMF in massive galaxies we require spectra with significantly 
higher S/N. Note however that such IMF sensitivity might 
be masked by similar index variations due to [Na/F e] abun - 
dance ratio, as shown in IConrov fc van Dokkuml l|201ll) . 
These authors studied the sensitivity of this feature to both 
[Na/Fe] and [a/Fe] abundances in SSP models of 13 Gyr. 

As it can be seen from Fig. [5] and Fig. [TS] this fea- 
ture strengthen with increasing IMF slope, whilst the Call 
triplet at ^860 A we akens. This extent has been n o ted b y 
ISchiavon et"afl (|200d ) and IConrov fc van Dokkuml l|201ll) . 
In the latter it is shown that an increase in [Na/Fe] causes a 
decrease in the strength of the Call triplet, as Na is a ma- 
jor contributor t o the electron pressure in late-type stars. In 
IVazdekis et~afl (|2003h we fully characterize the behaviour 
of the Call triplet in SSPs. Therefore these two features 
can be used to propose an index-index diagnostic diagram 
to constrain the IMF. Fig. [16] shows the resulting model 
grids for two different ages (2 and 12.6 Gyr) for models with 
Unimodal IMF. A range of IMF slopes and metallicities is 
shown. All the index measurements are performed on the 
LIS-14.0A system. These grids allow us to distinguish the 
IMF slope if the age and the [Na/Fe] abundance ratio are 
properly constrained. 
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Figure 14. Broadening correction for the newly defined 
NaI8200A index, ANaI8200A/NaI8200A, in SSPs with Unimodal 
IMF and varying slopes. The solid lines become thicker with in- 
creasing slope. The upper panel shows the results for a 12.6 Gyr 
model, whereas in the lower panel we show the results for a 
2 Gyr, the two for solar metallicity. We have broadened the 
model spectra by convolving with gaussians from the nominal 
resolution, i.e. cr=39 km s" 1 (FWHM=2.51 A) at 8200 A, up to 
cr=360 km s — 1 , in steps of 30 km s _1 . 



Figure 16. NaI8200A vs. CaT* diagnostic diagram to constrain 
the IMF slope. We adopt a U nimodal IMF for all the models. 
The CaT* index is defined in ICenarro et alj j200ll) . The solid 
lines indicate models with constant IMF slope, the thicker the line 
the higher the slope (jj. =0.3,0.8,1.0,1.3,1.5,1.8,2.0,2.3). The dot- 
ted lines represent models of constant metallicity ([Z/H]=-0.71,- 
0.40,0.0,0.22). The metallicity and IMF slope values are quoted 
in the left panel, where we show the grid corresponding to models 
of 12.6 Gyr. In the right panel we show the grid for 2 Gyr. All the 
index measurements are performed in the LIS-14.0A system. 



[Z/H] 

0.22 

0.00 

-0,40 




Age(Gyr) 

Figure 15. Behaviour of the NaI8200A index as a function 
of age for SSP models of different metallicitics, as indicated in 
the legend, for a Unimodal IMF with slope fj, = 1.3 (Salpctcr) 
and fi = 2.3, plotted in black and green lines, respectively. All 
the index measurements are performed in the LIS-14.0A system 
(FWHM=14A). 



4.3 Galaxy spectrophotometric surveys 

The use of multi-filter photometric surveys to determine 
SEDs and redshifts with hig h enough level of accuracy 
(SPSS, see also C OMBO- 17: IWolf et all |200S|; COS MOS: 
lllbert et al.1 120091 : ALHAMBRA: iMoles et~afl 12008ft has 
opened a new way to analyze the stellar populations of 
galaxies at different z-values and in different environments, 
thus allowing for the study of the evolution of galaxies and 
cosmology using huge samples. 

One of the major advantages of this kind of surveys is 
the fact that they provide low resolution spectroscopy for 
every pixel of the sky. With low resolution spectrophoto- 
metric data, spectral fitting techniques over the full spectral 
range are mandatory to exploit as much as possible the in- 
formation in the data. A unique characteristic of this type 
of data is the fact that the photo-spectra constructed on the 
basis of single narrow-band filter imaging does not suffer 
from flux calibration systematic uncertainties, unlike stan- 
dard spectroscopy. Every single point of the photo-spectrum 
-i.e. every filter- is observationally independent of the rest of 
the photo-spectrum, so the resulting SED is not affected by 
large scale systematics in the relative flux calibration (hence 
in the SED colors). With this unique advantage in mind, for 
the proper analysis of the SEDs of galaxies and stars it is cru- 
cial to use template stellar libraries with extremely accurate 
flux calibration, like MIUSCAT. In this sense, MIUSCAT, 
and all the SED SSP models derived from this library, are 
perfectly suited for the analysis and interpretation of op- 



14 



Vazdekis et al. 



M32 i-pas 




4000 5000 6000 7000 8000 9000 



X(k) 



comparisons. To facilitate the exploitation of these models 
we have integrated these new predictions within our recently 
developed website (http://miles.iac.es). This web provides 
the necessary support for using these model SEDs as well as 
those predictions based on MILES, CaT and Jones (1999) 
stellar library. 

All these models as well as the MILES and CaT stellar 
libraries can be retrieved and handled according the require- 
ments of the users. The webtools include the transformation 
of the spectra to match the instrumental set-up of the obser- 
vations (spectral resolution and sampling), measurements of 
line-strength indices and synthetic magnitudes derived from 
the spectrum. Further descriptions of all these applications 
are described in Paper II. 



Figure 17. Spectral fitting of M32 using the MIUSCAT SSP 
SEDs. The spectrum of the galaxy at the J-PAS resolution is 
plotted in red. The best MIUSCAT SSP fit (Age = 4.6 ± 1.1 Gyr; 
[Z/H]= 0.04 ± 0.08 dex) to the spectrum of M 32, as derived from 
a standard \ 2 minimization technique, is plotted in blue. The 
residuals are shown in the lower panel with the same scale. See 
the text for more details. 



tical spectrophotometric data, provided the accurate flux 
calibration of the library stars and the full optical spectral 
coverage. 

To illustrate a possible use of these new models, we 
focus on a test-case bas ed on the Javalambre-PAU Astro- 
phys ical Survey (J-PAS: iBenftez et al.l [20091 ; ICenarro et al.l 
2010), a new survey with 54 narrow-band (100 A) contigu- 
ous optical filters (AA ~3500 - 9500 A) under preparation 
at the Observatorio Astrofi'sico de Javalambre, which will 
provide low resolution (R~30) spectrophotometric data for 
hundred million galaxies. Figure [17| illustrates the best SED 
fitting derived for the integra ted spectrum of M 32 from 
iBica. Alloin. Sz Schmidt! (| 19901) , taken from the compilation 
of lSantos et all |2002l) . using the MIUSCAT SSP models as 
input templates. M 32 and the template spectra have been 
convolved with the J-PAS filters to simulate a real case. 
It is clear from the figure that the best fit, derived from 
a standard x' 2 minimization technique, reproduces well the 
observed spectrum at both low and high frequencies. The 
obtained residuals are shown in the lower panel. Note the 
telluric absorption still present in the data in the spectral 
range AA 7000-8000 A. The best solution corresponds to a 
MIUSCAT SSP model of 4.6 ± 1.1 Gyr and around solar 
metallicity (0.04 ± 0.08 dex), in good agreement with re- 
sults based_OTi_niuch_l2igh^ 

fe.g.lVazdekis fc Arimotolll999l :l Schiavon. Caldwell, fc Rose] 
2004). 



5 THE WEB TOOL 

The extension of the spectral range of our models based on 
the composite MIUSCAT stellar spectra allows us to cover 
the gap left by our previous model predictions based on the 
MILES and CaT libraries, as well as to extend blueward 
and redward these libraries the wavelength coverage. There- 
fore this extension represents a significant improvement that 
could be useful for many applications as well as for model 



6 SUMMARY AND CONCLUSIONS 

We have extended the spectral coverage of our stellar 
population synthesis models by combining the pred ictions 
based on the MILE S (|Sanchez-Blazquez et alj.l200tj) CaT 
jCenarro et al.ll200lh and Indo-U.S. (|Valdes et al.ll2004l ) em- 
pirical stellar spectral libraries. For this purpose we com- 
bined the stellar spectra of these three libraries for a sub- 
sample of 432 stars of the Indo-U.S. database, with no gaps 
in the relevant spectral regions and with no significant tel- 
luric absorption residuals. However the whole Indo-U.S. stel- 
lar database was taken into account for matching the stel- 
lar parameters of this library to the MILES/CaT system. 
Due to the limited flux-calibration quality of the spectra of 
the Indo-U.S. library, we used the MILES and CaT spec- 
tra as a reference. We scaled these two spectral ranges ac- 
cording to empirical relations between the stellar parameters 
and the Johnson- Cousin V — I colour, which are based on 
exten sive photomet ric stellar libraries (mainly I Alonso et al.l 
Il99tj for dwarfs and [A] onso et al. I l 19991 for giants). The same 
prescriptions were employed to co nvert the theoretic al pa- 
rameters of the stellar isochrones (jCirardi et al. 2000) that 
feed the models to observable fluxes and colours. For the se- 
lected Indo-U.S. stars that are lacking in either the MILES 
or CaT databases, or in the two, we synthesized the corre- 
sponding stellar spectra for the same atmospheric parame- 
ters using the MILE S and CaT interpolators d escribed in 
(|Vazdekis et aUl201Ch and (| Vazdekis et al.ll2003h . which are 
also employed to generate the required spectra to populate 
the isochrones during the SSP SED calculations. 

Appart of filling-in the spectral gap between the MILES 
and CaT libraries the Indo-U.S. stellar spectra allowed us 
to slightly extend bluewards MILES and redwards CaT the 
spectral ranges of these two libraries. Finally the spectral 
ranges covered by the Indo-U.S. and CaT stellar libraries 
have been smoothed and resampled to match the values of 
MILES. The resolution is therefore constant along the whole 
spectral range, i.e. 2.5 A (FWHM), with 0.9 A per pixel. The 
spectral coverage of the composite stellar spectral range is 
AA 3464.9 - 9468.8 A. 

The composite MIUSCAT stellar spectra were imple- 
mented in our models to compute SSP SEDs based on these 
spectra. These models were then combined with the pub- 
lished model SEDs for the M ILES ([Vazdekis et alJl201Ch and 
CaT ll Vazdekis et alj [20031 ) following a similar approach to 
that employed for the MIUSCAT stars, keeping the SSP 
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flux scaling for these two spectral regions. Therefore in the 
composite MIUSCAT SSPs the SEDs in the MILES and 
CaT spectral ranges are identical to the originally pub- 
lished models, with the exception of the resolution of the 
CaT models, which have been smoothed to match that of 
MILES. We estimated the reliability of the resulting model s 
according to the method described in IVazdekis et~afl |20ld) , 
which takes into account the atmospheric parameter cover- 
age of the 432 stars included in the MIUSCAT stellar library. 
The new model SEDs are reliable in the metallicity range 
(-0.71 < [Z/H] < 0.22) for all the ages (0.063 - 17.8 Gyr) 
and IMF shapes (Unimodal, Bi modal, Kroupa Univers al and 
Kroupa Revised, as described in lVazdekis et all (|2003r i) and 
IMF slopes (0.3^/i^2.3) that are considered in this work. 
Broad band colours are also acceptable outside these ranges. 

The new models allow us to measure colours involving 
the Johnson-Cousins and SDSS broad-band filters from U 
to / and from it to i, respectively. To measure the magni- 
tude in the U filter we compute a correction factor for the 
missing flux, as this filter extends bluewards the MIUSCAT 
spe ctral ranee. Th is correction is computed with the aid of 
the IPickled Jl998|) low-resolution stellar library. The good 
flux-calibration quality of the resulting model spectra al- 
lows us to measure accurate colours in the covered spectral 
range. The derived SSP colours are in very good agreement 
with the corresponding predictions obtained on the basis of 
empirical relations between colours and stellar parameters, 
which come from extensive photometric libraries. 

These model SEDs can be particularly useful for galaxy 
surveys, allowing to study the behaviour of the model SEDs 
with relevant stellar population parameters for a variety of 
filter responses. 

Unlike models based on theoretical stellar atmospheres 
our SSP models provide good fits to the colours of the glob- 
ular clusters of the Milky Way for the ages and metallicities 
that are derived from CMD studies. More extensive com- 
parisons with globular cluster data and galaxies are shown 
in the second paper of this series (Paper II). We show the 
dependence of the colours as a function of age, metallicity 
and IMF. We find that colours involving redder filters show 
greater sensitivities to the IMF slope. 

We also find various absorption line features and 
molecular bands that are very sensitive to the IMF slope 
throughout the whole spectral range covered by the MIUS- 
CAT model SEDs. Among the most prominent features 
we find that the Nal doublet at ~8200 A increases its 
strength with increasing IMF slope, whereas the neighboring 
Call triplet arou nd ~8600A, decreases, confirming previ- 
ous findings (e.g., Schiavon ct al. 2000i; IVazdekis et al.ll2003l ; 



IConrov fe van Dokkum 2011). Note that the Nal feature is 
not covered by the model SEDs based on the MILES and 
CaT stellar spectral libraries. 

With the aid of the MIUSCAT SSP SEDs we propose 
a new index definition for the Nal doublet at ~8200 A, 
NaI8200A, which has significant advantages over previous 
index definitions for this feature, particularly for low veloc- 
ity dispersion stellar systems. NaI8200A varies with the age 
and the metallicity and it is particularly sensitive to the IMF 
slope. We prop ose a NaI8200A - Ca T* diagram, the latter 
index defined in lCenarro et al.l (|200ll ). to distinguish the ef- 
fects of the IMF slope if the age and the [Na/Fe] abundance 
are properly constrained. 



Finally we show the utility of these new models for 
survey-oriented spectrophotometric applications. Because of 
the wide spectral coverage and accurate flux calibration of 
the MIUSCAT library stars, the SSP models are perfectly 
suited for being used as templates for exploiting galactic 
spectrophotometric data, deriving reliable estimations of 
ages and metallicities even for low resolution data. 

The MIUSCAT model SEDs, as well as all our model 
predictions, can be downloaded from the miles website 
(http://miles.iac.es), where web-based facilities are avail- 
able for downloading and handling these SEDs. These user- 
friendly tools also allow to obtain the MIUSCAT SEDs, with 
the corresponding line-strength indices and colours, for a 
variety of SFH parametrizations as well as for user-defined 
SFHs. 
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